Prospects for a mHz-linewidth laser 
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We propose a new light source based on having alkaline-earth atoms in an optical lattice col- 
lectively emit photons on an ultra-narrow clock transition into the mode of a high Q-resonator. 
The resultant optical radiation has an extremely narrow linewidth in the mHz range, even smaller 
than that of the clock transition itself due to collective effects. A power level of order 10 _12 VK is 
possible, sufficient for phase-locking a slave optical local oscillator. Realizing this light source has 
the potential to improve the stability of the best clocks by two orders of magnitude. 
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Time and frequencies are the quantities that we can 
measure with the highest accuracy by far. From this 
fact derives the importance of clocks and frequency stan- 
dards for many applications in technology and fundamen- 
tal science. Some applications directly relying on atomic 
clocks are GPS, synchronization of data and communica- 
tion networks, precise measurements of the gravitational 
potential of the earth, radio astronomy, tests of theories 
of gravity, and tests of the fundamental laws of physics. 

With the advent of octave spanning optical frequency 
combs [l|, Q it has become feasible to use atomic transi- 
tions in the optical domain to build atomic clocks. Opti- 
cal clocks based on ions [3|and ultracold neutral atoms 
confined in optical lattices [4j have recently demonstrated 
a precision of about 1 part in 10 15 at one second and a 
total fractional uncer tainty of 1CT 16 [|] or below [9] , sur- 
passing the primary cesium microwave standards [a, 0] . 

The state-of-the-art optical atomic clocks do not 
achieve the full stability that is in principle afforded by 
the atomic transitions on which they are founded. These 
transitions could have natural line-Qs of order 10 18 , ex- 
ceeding the fractional stability of the clocks by a factor 
of ~ 100. The main obstacle that prevents us from reap- 
ing the full benefit of the ultra-narrow clock transitions 
is the linewidth of the lasers used to interrogate these 
transitions. These lasers are stabilized against carefully 
designed passive high-Q cavities and achieve linewidths 
< 1 Hz, making them the most stable coherent sources of 
radiation. It is mainly the thermal noise of the reference 
cavity mirrors that prevent a further linewidth reduction 
and substantially reducing this noise is hard 

An elegant solution to these problems would be to di- 
rectly extract light emitted from the ultra-narrow clock 
transition 0. That light could then be used as an opti- 
cal phase reference, circumventing the need for an ultra 
stable reference cavity. Unfortunately, the fluorescence 
light emitted on a clock transition is too weak for practi- 
cal applications. For instance, for 10 6 fully inverted 87 Sr 
atoms the power of the spontaneously emitted light is of 
the order of 10" 16 W. 

The key observation that motivates this work is that 
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FIG. 1: (Color online) (a) Experimental scheme and (b) the 
level structure. The ultra-narrow 3 Po~ 1 So transition indi- 
cated by the red arrow is the laser transition. The repumping 
lasers are indicated by the blue arrows, spontaneous decay 
into the 3 P states by the gray arrow, and the Raman side- 
band cooling lasers by the black arrows. The n quantum 
numbers indicate the vibrational levels of 3 Po and 3 P2- 



if we could coerce the ensemble of atoms to emit the 
energy stored in them collectively rather than individ- 
ually, the resulting power of order 10 -12 W would be 
large enough to be technologically relevant. We show 
that such collective emission of photons can indeed be ac- 
complished if the atoms are located inside a high-Q cav- 
ity. The collective interaction between atoms and cavity 
fields HHH [3 El as well as lasers based on emission 



from microscopic atom samples 1J, Il5l . Il6l . |17| are of 
great theoretical and experimental interest in quantum 
optics. 

Such a laser would operate in a strikingly different 
regime from conventional lasers: the cavity relaxation 
rate exceeds the atomic relaxation rates by many orders 
of magnitude. This system is thus an extreme case of a 
bad-cavity laser. Such bad cavity lasers have been stud- 
ied in the past 3, H, 2(| , although in those papers the 
separation of timescales was not nearly as large as is con- 
sidered here. In the microwave domain masers are prime 
examples of a bad cavity laser. Also in the previous sys- 
tems the cavity field contains macroscopic numbers of 
photons while in our system the field occupation number 
can be < 1. The extremity of the current system makes 
it necessary that we revisit the foundations of the the- 
ory to obtain reliable predictions. We will focus on the 
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example of 87 Sr confined in a lattice potential to make 
the presentation clearer and more concise. Our results 
however are general and can be easily translated to other 
alkaline-earth-like atomic systems with similar structure. 

The general setup is shown schematically in Fig. [T] We 
consider N ultra-cold two level atoms with transition fre- 
quency Lu a confined in an optical lattice, e and g corre- 
spond to the two clock levels, 3 Po and 1 5'o in 87 Sr. The 
atoms have a spontaneous decay rate 7 and we model all 
inhomogeneous processes by an effective relaxation rate 
T 2 _1 for the atomic dipole. The laser transition is cou- 
pled to a single mode of a high Q-optical resonator with 
resonance frequency lu c and linewidth k. We assume that 
the atoms are held in place by the external optical lat- 
tice potential exactly in such a way that all atoms cou- 
ple maximally and with the same phase to the specific 
cavity mode. Repumping lasers resonantly drive transi- 
tions from the 1 S'o to the 3 Pi transition and from there 
to the 3 Si state. The repump lasers optically pump the 
atoms into the 3 Pq and 3 p2 states. A Raman transition 
from 3 P to 3 P 2 via 3 Si is used to implement side-band 
cooling to the vibrational ground state and to optically 
pump all atoms to the 3 Po state thus providing inversion 
for the laser transition. We summarize all these repump 
steps by an effective repump rate w. We have estimated 
AC Stark shifts incurred by these repump and cooling 
lasers and we find that for repump rates up to around 
w ~ 10 3 s _1 the AC stark shift induced fluctuations of 
the atomic transition frequency is negligible at the mHz 
level, if the laser intensity is stabilized to 1%. With the 
continuous cooling in place we are justified in assuming 
that the atoms are in the ground state of the lattice. 

A discussion of the relative and absolute magnitude 
of the different rate constants is in order. For 87 Sr, the 
linewidth of the doubly forbidden intercombination line 
is 7 w 0.01 s _1 . The inhomogeneous lifetime in the most 
recent generation of lattice clock experiments has been 
pushed to T2 ~ 1 s [2l|. The effective repump rate w 
can be widely tuned from to values of order 10 4 s _1 
and beyond. The cavity parameters enter the physics 
through the cooperativity parameter C — fl 2 / (ivy) which 
is independent of the length of the cavity. For definitc- 
ness we consider a cavity with an effective mode volume 
V e s = (lmm)7rx (50/mi) 2 . Because of the extremely weak 
dipole matrix element of the inter combination transition 
of order 10 _5 eao, with e the electron charge and ao the 
Bohr radius, this leads to a vacuum Rabi frequency of 
order Q ~ 37s . In this setup the cavity decay rate is 
by far the largest time scale. For a finesse of T = 10 6 
and the above cavity parameters the cavity linewidth is 
k = 9.4 x 10 5 s -1 , and thus C « 0.15. 

The coupled atom-cavity system can be described by 
the Hamiltonian 
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In this formula, fl — hr 1 y/fvjj c / (2eoVes), eo the vacuum 
permittivity, and a and at are bosonic annihilation and 
creation operators for photons in the laser mode. We 
have introduced Pauli matrices ctJ = |ej)(ej| — \gj)(gj\ 
and a J = (ojT)t = \9j)(^j\ f° r the jth atom. 

We take the various decay processes into ac- 
count by means of the usual Born-Markov master 
equation for the reduced atom-cavity density ma- 
trix p, dp/dt = (ift) _1 [P, p] + C[p], with the Li- 
ouvillian C[p] — £ C avit y [/6] + £ spont. [p] tinhorn. [Pi 
^repump [p] ■ The Liouvillian for cavity decay is 
£cavity[/5] = — n/2(d^dp + pata — 2apat), the sponta- 
neous decay of the atoms is described by £ S pont[p] = 

-7/2 D^li {®t°I P + P°t°j ~ P°t) ' and tlie Liou " 
villian for the inhomogeneous life time is A n hom.[p] = 

l/(2T 2 )Y,f =1 (vjPVj ~ P)- The Liouvillian for the re- 
pumping is identical to that for the spontaneous decay 
with the replacements 7 — > w, &J — > o - ^, and a~j — * <jJ . 

An important aspect of this system that is born out 
by the master equation is that the coupling of the atoms 
to the light field is completely collective. The emission of 
photons into the cavity acts to correlate the atoms with 
each other similar to the case of ideal small sample super- 
radiance 23, HI, leading to the build-up of a collective 
dipole. The locking of the phases of the dipoles of differ- 
ent atoms gives rise to a macroscopic dipole that radiates 
more strongly than independent atoms. The macroscopic 
dipole is also more robust against noise from decay pro- 
cesses and repumping, leading to a reduced linewidth. 

We have verified that the system settles to steady state 
much faster than the anticipated total operation time, 
provided that the repump rate is not too close to the 
laser threshold derived below. For the representative ex- 
ample parameters used below, the relaxation oscillations 
decay after a time < Is while the total operation time 
will typically be > 1 minute. We therefore focus en- 
tirely on the steady state behavior in this Letter. To 
find the steady state we introduce a cumulant expan- 
sion to second order for the expectation values of system 
observables (24|. We denote raw expectation values by 
(. . .} and cumulant expectation values by (. . .) c . In our 
model all expectation values, cumulant or raw, are com- 
pletely symmetric with respect to exchange of particles, 
for instance (<5j) c = (cti) c and {d^aj) c — (d^a^) c for 

all j, (c r f(Tj)c = (vi&2)c f° r au i h e tc- In this 
formalism where we explicitly keep track of higher order 
correlations, the total phase invariance of the system is 
not broken and we have (a) — (at) — (af) = 0. The 
only non-zero cumulant of the first order is then the in- 
version (<7i) c = (o"i). This cumulant couples to second 
order cumulants through the atom-field coupling, 
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where do = (w — "f)/(w + 7). The atom-field coherence 
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evolves according to 
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where 5 = lu c — u> a . We have estimated the third order 
cumulant in square brackets to be much smaller than the 
other terms and therefore neglect it. The second and 
third terms represent the exchange of energy between 
cavity field and a single atom. They are non-collective 
in nature. The last term describes the coupling of the 
atom-field coherence to the collective spin-spin correla- 
tions which locks the relative phase between atoms and 
field to the phase of the macroscopic atomic dipolc. 
The spin-spin correlations evolve according to 



dt 



= -(w + 1 + 2T 2 - 1 )(a+^) c 



(4) 
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where we have dropped the small third order cumulants 
of the type (a^of d-,7} c . To close the set of equations we 
also need the equation for the mean photon number, 



d{a'a) c 
dt~ 



= -K(a t a) c + -— ((a t a 1 ) c - (&fa) c ). (5) 
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We consider the steady state of this system for the case 
6 = by setting the time derivatives in Eqns. (HJ)-© to 
zero. The resulting algebraic equations can be solved 
exactly. Simple approximate results can be obtained in 
certain limits on which we will focus in our discussion 
to explain the underlying physics. The numerical results 
reproduced in the figures are based on the exact solutions 
and agree well with the approximate treatment. 

Our first goal is to understand the role of the collective 
effects. Neglecting all decay constants but k in Eq. ([3]), 
keeping only the collective term proportional to (o"i~ct^~) c , 
and approximating N — 1 w N we find for the atom-field 
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serting this result into the steady state equation for the 
inversion determines the saturated inversion and plug- 
ging that and (a^a^) into the steady state equation for 
the spin-spin correlations yields the central equation 



= (<r+a 2 -) c -r + doN-fC - 2 



w + 7 



(at &2 )c 



(6) 

where T = 7 + W + 2/T2 is the total relaxation rate of 
the atomic dipole. The solution for (o" 1 "6'^) c correspond- 
ing to the term in parenthesis is the physically stable 
solution. The laser threshold is the pump rate at which 
the gain dgN^/C overcomes the losses T. In the limit 
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FIG. 2: (Color online) Power as a function of pump rate w and 
atom number N. The rapid build up of power above threshold 
w ~ 7 can be seen as well as the decrease of emitted power 
for too strong a pump. The dashed line shows the boundary 
of the region of collective emission as determined by the zero 
of the term in parenthesis in Eq. Here, 7 = 0.01 s , 

1/T 2 = Is" 1 , n = 37s" 1 , and k = 9.4 ■ 10 5 s" 1 . 



T/^NC) — > this condition turns into w > 7. At 
threshold the pump overcomes the atomic losses which 
is in contrast to conventional lasers where threshold is 
obtained when the pump overcomes the cavity losses. 

At threshold the spin-spin correlations change sign sig- 
nifying the onset of collective behavior. Interestingly 
the spin-spin correlations change sign again at a larger 
pump rate above which the atoms return to normal non- 
collective emission. This upper threshold comes about 
because do eventually saturates at 1 while the pump in- 
duced noise grows with w. Setting do = 1 and neglecting 
all atomic noise sources other than w we find for the 
maximum pump rate 



NC-f 



(7) 



Above this threshold the pump noise destroys the coher- 
ences between different spins faster than the collective 
interaction induced by the light field can establish them. 

A minimum number of particles is necessary to have 
collective behavior. Below this critical number, (a^ <J^)c 
is never positive and hence the spin-spin correlations are 
never collective. The critical particle number can be 
estimated from Eq. (|6|). Assuming T27 *C 1, we find 
iVcrit = 2/(^7^2)- Physically, this equation means that 
there must be enough particles for the system to be in 
the collective strong coupling regime. 

The collective versus non-collective behavior is nicely 
illustrated by the outcoupled laser power as a function of 
w and N (Fig. [2]). Above threshold the outcoupled power 
rapidly increases until the collective dipole is destroyed 
at the second threshold w max - 

Equation ([6|) also allows us to determine the maximum 
spin-spin correlation of (o'J l "o' 2 _ ) c = 1/8 which is obtained 
for the pump rate w opt = NC^ /2. At this pump rate the 
laser power reaches its maximum of 



hw a N 2 Cj/8. 



(8) 
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The scaling of that power with the square of the num- 
ber of atoms underlines the collective nature of the emis- 
sion. Remarkably, this power is only a factor 2 smaller 
than the power expected for perfect super-radiant emis- 
sion from the maximally collective Dicke state at zero 
inversion. Figure [5] shows that an outcoupled power of 
order 10 -12 W is possible with 10 6 atoms. 

From the perspective of potential applications the most 
striking feature of this laser is its ultra-narrow linewidth. 
To find the spectrum we use the quantum regression the- 
orem to find the equations of motion for the two time 
correlation function of the light field (aJ(t)a(0)}. This 
correlation function is coupled to the atom-field correla- 
tion function (<7+(t)a(0)). Factorizing (a z (t)a^ (t)a(O)) « 
(<7 Z (t)) (a* (t)a(0)) we arrive at the closed set of equations 



di 



<at(t)a(0)> 
(* + (t)a(0)) 



-k/2 



2 

-r/2 



(at(i)a(O)) 
(*+(t)o(0)) 



(9) 

The initial conditions are the steady state solutions dis- 
cussed above. Laplace transforming the solutions yields 
the spectrum which is Lorentzian with linewidth Av for 

UJ a = LU C . 

For our example parameters Av is shown in Fig.[3j The 
leftmost dashed line in that figure is 7 corresponding to 
the threshold for collective behavior. When the pump 
strength w passes through that threshold the linewidth 
gets rapidly smaller with increasing w. When w reaches 
I/T2, indicated by the second dashed line, essentially all 
atoms are phase locked together. From that point on the 
pump noise due to w grows in proportion to the size of 
the collective spin vector. Therefore the linewidth is ap- 
proximately constant. Making similar approximations as 
in the steady state calculations we obtain the estimate for 
the minimum laser linewidth, Av = C7. That estimate 
agrees well with our numerical results. It is important 
to note that the parameters for achieving the maximum 
outcoupled power Eq. 1(5)) and the minimum linewidth 
are compatible with each other. For our example param- 
eters the estimate yields a linewidth smaller than the 
homogeneous linewidth of the atomic clock transition. 
We have studied the dependence of the laser frequency 
on the atom-cavity detuning and we find that this nar- 
row linewidth can be observed if the cavity resonance 
frequency is stable at the 1 kHz level which is relatively 
easy to achieve experimentally. When w increases beyond 
io maxl indicated by the third dashed line, the collective 
dipole is destroyed and the linewidth increases rapidly 
until it is eventually given by w. 

Future research is targeted at fully understanding the 
recoil effects and the detailed nature of the joint atomic 
and field state. 

We gratefully acknowledge stimulating discussions 



with D. Kleppner, J. K. Thompson, and J. Cooper. This 
work was supported by DFG, DARPA, NIST, DOE, and 
NSF. 

o 5 

o 3 „ 






-7— 

















































1(T 3 KT 1 10 1 10 3 10 : 



FIG. 3: (Color online) (a) Linewidth vs. w and N. The white 
dashed lines indicate (from left to right) the spontaneous de- 
cay rate 7, the inhomogeneous relaxation rate I/T2, and the 
maximum pump rate u) max (Eq. 10). Parameters as in Fig. [2] 
(b) is a cut through figure (a) for N = 10 6 atoms. 
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